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Abstract 

We study the response of the transmission eigenvalue spectrum of disordered 
metallic conductors to an arbitrary external perturbation. For systems without 
time-reversal symmetry we find an exact non-perturbative solution for the two-point 
correlation function, which exhibits a new kind of universal behavior characteristic 
of disordered conductors. Systems with orthogonal and symplectic symmetries are 
studied in the hydrodynamic regime. 



PACS numbers: 72.10.Bg, 05.60,+w, 05.40.+j 



1 



Since its first application to the statistics of level spacings in the energy spectra 
of complex nuclei [1] , the random matrix theory of Wigner, Dyson, Mehta and others 
has been developed into a powerful non-perturbative method to describe a number 
of local statistical properties of many complex physical systems, such as disordered 
metallic grains [2] , quantum chaotic systems [3] , random surfaces [4] and disordered 
conductors [5]. Random matrix ensembles can be conveniently derived, for any 
given spectral density, by means of a maximum-entropy requirement, which in turn 
implies that the local spectral fluctuations are determined solely by some universal 
symmetries of the system characterized by the parameter f3: f3 = 1 for systems with 
time-reversal symmetry and without spin-orbit scattering, j3 = 2 for systems without 
time-reversal symmetry and j3 = 4 for systems with time-reversal symmetry in the 
presence of spin-orbit coupling. It is precisely this universal behavior that accounts 
for the success of the many applications of random matrix theory. 

It is often difficult to find a microscopic justification for the maximum-entropy 
ansatz that generates the random matrix ensembles. For the particular case of dis- 
ordered metallic grains, however, Efetov [6] has derived a microscopic theory which, 
in the limit where a particle can diffuse throughout the sample, yields precisely the 
same two-point correlation function for pairs of energy levels obtained by random- 
matrix theory. 

Very recently, Altshuler, Szafer and Simons [7] have extended Efetov's theory 
to study spectrum correlations of disordered metallic grains in the presence of an 
arbitrary external perturbation. The most striking conclusion of their work is that 
the two-point correlation function for pairs of energy levels at different values of the 
external perturbation becomes universal if the energy levels and the perturbation 
parameter are appropriately rescaled. In the hydrodynamic regime the "velocity- 
velocity" correlator C(Su) [7] has been shown to have a remarkably simple form 

cm = m 

where Su is the perturbation parameter. Going beyond the hydrodynamic regime, Si- 
mons and Altshuler [7] have derived an exact non-perturbative expression for C(5u), 
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which after appropriate rescaling of the physical parameters becomes a universal 
function which is believed to be a new quantum mechanical characterization of sys- 
tems that obey the Wigner-Dyson statistics [8] of consecutive level spacings. 

Subsequently, Simons, Lee and Altshuler [9] have uncovered a number of correla- 
tions between these results and other problems. In particular, they have shown that 
C(5u) can also be derived from the periodic Sutherland Hamiltonian [10], which 
is known to be equivalent to the Dyson Brownian motion model [10]. The case 
(3 = 2 is particularly simple, since the interaction term in the Sutherland model 
vanishes and an exact solution for the two-point correlation function has been cal- 
culated. A detailed account of the application of Dyson Brownian motion model to 
describe parametric correlations in the spectra of disordered metals has been given 
by Beenakker and Rejaei [11]. 

The universality discussed by the above authors occurs when the eigenvalues of 
the random matrix are rescaled at the bulk of the spectrum. However, large random 
matrices are known to have eigenvalue densities limited to a finite support and it 
has been found [12,13] that, in the absence of external perturbation, different kinds 
of universal behavior emerge depending whether the eigenvalues are rescaled at the 
endpoints of the support or at the bulk of the spectrum. It is therefore natural to 
expect a similar scenario when an adiabatic external perturbation on the eigenvalues 
is turned on. The purpose of this letter is to demonstrate that this is indeed what 
happens. In particular, we show that the parameter dependent two-point correlation 
function for pairs of eigenvalues belonging to the Laguerre ensemble gives rise to a 
new universal behavior when rescaling takes place at the hard edge of the support. 
Our main result (Eq. (4)) is a direct consequence of the Bessel kernel [13] in the 
same way as the result of Refs. 7, 9 and 11, for f3 = 2, is a direct consequence of 
the sine kernel. Since the Bessel kernel is a local property which is known to appear 
at the hard edge of the support of the eigenvalue density of many random matrix 
ensembles, we believe that Eq. (4) is of equal general validity. 

The Laguerre ensemble has been applied to problems in several areas of theoret- 
ical physics ranging from nuclear physics [14], QCD [13] and disordered conductors 
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[5]. In this letter we shall be mostly concerned with the random matrix theory 
of disordered conductors. The fundamental difference between the random matrix 
ensembles appropriate to the description of disordered conductors and the more fa- 
miliar gaussian ensembles, used for modelling disordered metallic grains, is that one 
is not concerned with the energy levels, — oo < Ei < oo, of a random Hamilto- 
nian, but rather with the parameters Aj = 1/Tj — 1, where Tj are the eigenvalues of 
ttf (t is the transmission matrix), which because of flux conservation must satisfy 
< Aj < oo. The positivity constraint on the Aj's implies that the two-point corre- 
lation function is not translation invariant. Consequently, the variance of a general 
linear statistic [15] differs from that for gaussian ensembles, although the universal 
1//3 dependence is the same in both cases. 

It is interesting to compare random matrix theory with the microscopic theory 
of disordered mesoscopic conductors that has been developed by Lee and Stone and 
by Altshuler [16]. In this theory one is interested in mesoscopic fluctuations of the 
conductance, G, with the Fermi energy , Ep, or an applied magnetic field , B. This 
information is contained in the correlator F(5u) = (G(5u)G(0)) — (G(0)) 2 , where 
(. . .) denotes ensemble average and 5u is a small perturbation on E F or B. F(Su) has 
been shown to decay as a power law as a function of 5u/£ c , where S c is a correlation 
parameter which sets the typical scale of the spacings between the peaks and valleys 
in the values of G, for a given realization, as a function of 5u. This slow power law 
decay of F(5u) indicates that the system sustains some sort of long-range memory in 
the fluctuations of the conductance. This kind of behavior must be contrasted with 
the exponential decay of Poisson processes and can be regarded as a manifestation 
of "spectral rigidity", a well known phenomenon in the theory of random matrices. 

For the sake of mathematical tractability we shall consider the global maximum- 
entropy approach [5]. This approach has been shown [17] to be an exact description 
of the local statistical properties of levels that are close to zero on scales comparable 
with the mean level spacing. For levels away from the origin this approach provides 
an adiabatic approximation for the actual distribution which turns out to be an 
excelent approximation in the hydrodynamic limit. In the metallic regime [5] one 
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has the following probability distribution for the variables Xi xs \f\l 

Po '({Xi}) =C%U K - °?jf U x i a+1 ex P(-f CX ^ ( 2 ) 

i<j i=l 

where c is a system dependent constant and a is a parameter that as we shall 
demonstrate below labels different universality classes for the two-point correlation 
function (a = for disordered conductors). The most general case, where a > — 1, 
was first derived by Bronk [14] and characterizes the Laguerre ensemble of random 
matrices. 

We now state our results. Assuming that spectrum correlations are dominated 
by level repulsion, we propose a Brownian motion model (to be described below) for 
the evolution, with an external perturbation parameter 8u, of the eigenvalues {x^, 
which is consistent with the asymptotic "equilibrium" distribution Pq({xi}). For 
(3 = 2 we show that the resulting Fokker-Planck equation can be mapped exactly 
onto a free-fermion Hamiltonian. We calculate exactly the two-point correlation 
function defined as 

S(x, x', 8u) = $( x — Xi(u))5(x' — Xj(u + Su))^ — p(x)p{x'), (3) 

ij ^ 

where (. . ) eq stands for an arbitrary average calculated with the equilibrium distri- 
bution P^({xi}) and p{x) = (^2 i 8(x — Xi(u))) eq . After rescaling the eigenvalues at 
the hard edge of the spectrum we find 

rl poo 

S(x, x 1 1 8u) =7r 4 p(j ds ds / ss'xx'J ct (rcpoxs)J a (Trpox's)J a (rcpoxs') 

J a ( W V) exp(^^ 2 ( S 2 - S /2 )), (4) 

in which J a (x) are the Bessel functions, po is related to the average density p(x) 
through p(x) = 7r 2 Po Jo dssxJ^npQXs) and 7 _1 = {13/2) (xf(u)), where x\ = dxi/du, 
is the friction coefficient of the Brownian motion model. Note that our expressions 
for p(x) and S(x, x', 0) agree with known results [13]. The correlation parameter S c 
discussed in the introduction must have, by definition, the same dimension as 8u. 



Since Xi and po are dimensionless it is natural to introduce the definition £~ 2 = 
f3p {x 2 (u)} /2. A straightforward evaluation yields £ c = a/tTpo- 

Consider now the particular case of disordered conductors (a = 0). Following 
Ref. 7 we make the rescalings x = pox, x' = pox' , Su 2 = Su 2 P0/7 an d ■> $ u ) = 

S^rr, a;', 8u)/pq. We can see that «j>(x, <5w) is a function that is independent of any 
physical parameter, such as the Fermi velocity, the mean free path or the length of 
the sample, and therefore can be regarded as a universal characterization of such 
systems. Since we started from the global maximum-entropy ansatz, Eq. (4) is 
exact for disordered conductors only for levels that are close to the origin on scales 
of the order of the mean level spacing. For the more general case of the Laguerre 
ensemble, with a > —1, we can see that the parameter a labels new universality 
classes characteristic of systems with a hard edge in the spectrum. 

The "velocity- velocity" correlator defined as 

C(x, x', 8u) = Xi(u)xj (u + 5u)5(x — Xi{u))5{x — Xj(u + 5u))^ , (5) 

ij ^ 

is related to (4) by means of the identity d x d x >C(x 7 x', u) = —d 2 S(x 7 x', u). We shall 
consider {xi) eq = without loss of generality. 

In the hydrodynamic regime fluctuations in the level density are small [18]. 
Using the fact that p(x) w po for x 3> 1/po, we get C(x, x, 8u) w — 2/(f3n 2 8u 2 ), 
which is valid for arbitrary (3 and Su ^> u c . This is precisely the universal correlator, 
C(5u) of Eq. (1), found in Ref. 7. We can estimate u c if we recall that the integral 
of C(x, x, 8u) over all 8u vanishes by definition, and use the definition of 7 to get 
'J~ 1 PqU c ~ — f u C(x, x, Su)dSu. Evaluating the integral we find u 2 ~ 7/po- For 
8u < u c fluctuations in the level density are not small and C(x, x, Su) becomes, in 
general, an ct-dependent function. 

We shall now only outline our calculations. The mathematical details will be 
presented elsewhere [19]. Following Dyson's idea [10], we model the parametric de- 
pendence of the distribution P({xi}, r) of the eigenvalues {x^ by Brownian motion 
of iV classical particles at positions {xi(r)} moving in a viscous fluid with friction 
coefficient 7 _1 at temperature f3~ l . The parameter r is a fictitious "time", which 
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is related to the perturbation parameter Su by r = 5u 2 [9,11]. We require the ini- 
tial condition P({xi}, 0) = [| i S(xi — where are the transmission eigenvalues 
of a perturbation-free system. We also require that as r — > oo the system tends 
to a maximum-entropy configuration, therefore P({xi},oo) = Pq*({:e£}). The sim- 
plest description of a Brownian motion satisfying these requirements is given by the 
Fokker-Planck equation 

where 

i i i<j 

More general descriptions can be obtained by allowing for other sources of correla- 
tions in addition to the pure geometric effect of level repulsion (see Ref. 11 for a 
short discussion). As in the case of the conventional Dyson Brownian motion model 
we shall ignore these additional effects and assume that correlations are dominated 
by level repulsion. 

We proceed by writing (6) and (2) in terms of the variables Aj = x\ and then 
we perform the transformation P({AJ,t) = P Q ({A i }) 1 / 2 *({A l }, r). We find that 
^({Ai}, r) satisfies the Schrodinger equation in imaginary time —d^/dr = (H—B)^, 
where 



,, v- a a fny, , a\ , gtg-2) ^ 



(7) 



B = Nc(2 + P(N + a-l))/(2^), u = 4c/(^y) and n = cA,. Observe that for (3 = 2, 
TC becomes a non-interacting Hamiltonian, whose single electron eigenfunctions are 
the Laguerre functions <E>" (A). Using its asymptotics and taking N — > oo and c — > 0, 
such that 2\fWc — > 7T/9 , we can evaluate a;', 5m) exactly to find Eq. (4). A 
more detailed discussion of «S(a;, x', 5m) can be found in Ref. 19. We now discuss the 
hydrodynamic regime. 

The hydrodynamic limit can be calculated directly from (6) for arbitrary f3. 
Multiplying (6) by ^2 i S(x — Xj) and integrating over all x^s we get the following 



evolution equation for the average density 



dp(x,r) d 
7—!— PS 

dr dx 



d ( r°° 

p ^ x ' T ^~dx\J d yp(y^ T ) ln \ x2 - y 2 



C 2 
-X 

2 



(8) 



in which only the dominant contributions for large N were considered. The equi- 
librium density is p eq (x) ~ 2n~ 1 \ / Nc — > po = const, and thus tends to a constant 
once the limiting procedure introduced above is implemented. 

We are concerned with the diffusive regime, where the average dimensionless 
conductance, (G) eq: satisfies (G) eq 3> var(G) ~ 0(1). Since according to the Lan- 
dauer formula the average conductance is simply (G) eq ~ po, one expects fluctuations 
in Peq(x) to be small in this regime. Defining Sp(x, r) = p(x, r) — po, inserting into 
(8), neglecting contributions of second order in dp and taking the Fourier cosine 
transform of the resulting equation, we find Sp(k,Su 2 ) = Sp{k, 0) exp(— knSu 2 /£ 2 ). 
From (3) we get for the Fourier transform of S(x,x',Su) the result S(k,Su) = 
k/(2f3)exp(—kT:8u 2 /£c) : in which we have used Beenakker's result [15] for S(k,0). 
This equation is valid in the range k po for all Su. 

We are now in a position to calculate C(x,x,Su). Using the relation between 
C(x, x', u) and S(x, x', u) and the above expression for S(k, Su) we find C(x, x, Su) ~ 
—2/(f3n 2 Su 2 ), for Su 3> u c and x 3> l/po ; this result is just Eq. (1). 

Finally, we shall use our result for S(x, x', Su) to calculate parametric correla- 
tions of the conductance. Defining SG = G — (G) eq and F(Su) = (SG(Su)SG(0)) eq 
we find 

7T SU* 

Observe that when Su = 0, Eq (9) gives var(G) = 1/(8/3), which is the value found 
in Ref. 15. For Su < £ c one gets from (9) F(Su) ps 1/(8/3) (1-ttSu 2 /£ 2 ), a quadratic 
decrease, as found in Ref. 16. For Su 3> £ c we find F(Su) ps (Su/£ c )~ 4 , which is 
a power law decay as one would have expected from spectral rigidity and is also 
in qualitative agreement with Ref. 16. The exact value of the exponent found in 
Ref. 16 shows a nontrivial dependence on both the dimensionality and the nature of 
the perturbation, which cannot be obtained from the present random matrix model. 



o N -2 



F (Su) = ^[l + -^) . (9) 



The overall shape of F(5u) in Eq. (9), therefore, is in reasonable good qualitative 
agreement with microscopic calculations. 

In conclusion we have presented a consistent theory of parametric correlations 
in the Laguerre ensemble. As a particular case we have considered the problem of 
the response to an arbitrary external perturbation of the transmission eigenvalue 
spectrum of disordered metallic conductors. We believe that some of our interme- 
diate results, like Eqs. (6) and (7), might be of general interest in connection with 
applications of the Laguerre ensemble in other branches of physics. 
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